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ABSTRACT 

Starting from a sample of SDSS quasars appearing also in the 2MASS survey, we study 
the continuum properties of ^ 1000 objects observed in 8 bands, from NIR to UV. We 
construct the mean spectral energy distribution (SED) and compare and contrast the 
continua of radio loud (RLQ) and radio quiet (RQQ) objects. The SEDs of the two 
populations are significantly different, in the sense that RLQs are redder with power 
law spectral indices (q;r,lq) = —0.55 ± 0.04 and (aRgq) = —0.31 ± 0.01 in the spectral 
range between lO^**'^ and 10^^ '^''^ Hz. This difference is discussed in terms of different 
extinctions, different disc temperatures, or slopes of the non-thermal component. 

Key words: galaxies: active - galaxies: nuclei - quasars: general. 



1 INTRODUCTION 

A substantial effort has been dedicated to construct the 
spectral energy distributions (SEDs) for sizeable samples 
of quasars over the whole accessible range of the electro- 
magnetic spectrum (see e.g. Sanders et al. 1989; Francis et 
al. 1991; Elvis et al. 1994; Richards et al. 2006). However, 
only a few papers focus on the comparison between the SEDs 
of radio loud (RLQ) and radio quiet (RQQ) quasars. Elvis 
et al. 1994 propose an overall spectrum from a sample of 47 
quasars, divided in RLQs and RQQs, which shows that no 
distinction between the SEDs of the two subsamples is ap- 
parent in the range 100^ — lOOOA. As noted by these authors, 
the considered sample is biased towards X-ray and optically 
bright quasars. Some indication of a possible difference be- 
tween the NIR to optical colors of RLQs and RQQs in the 
2MASS catalogue was reported by Barkhouse & Hall (2001). 
Francis, Whiting & Webster (2000) found that the optical- 
NIR continuum is significantly redder in radio selected RLQs 
from the PKS Half- Jansky Flat-Spectrum survey than in op- 
tically selected RQQs from the Large Bright Quasar Survey 
(LBQS). These and other works (e.g. Kotilairien et al. 2007) 
indicate that radio loud objects are possibly redder than 
their radio quiet counterparts, but the samples seriously 
suffer of selection effects against red radio quiet quasar, be- 
cause RQQs are mostly collected from optical selected sam- 
ples (e.g. the LBQS, the first selection criterion of which is 
"blue color" of candidates; see Hooper, Impey & Foltz 1997). 
White et al. suggest that redder quasars from the Sloan Dig- 
ital Sky Survey (SDSS) are likely to be more radio-powerful 
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than bluer objects, and this is an indication that the red 
color of radio loud objects is not completely due to the bias 
against red RQQs in optically selected samples, since the 
two classes of objects derive from the same survey. 

Because of the importance of the distinction between 
RLQs and RQQs in the Unified Models of AGN, these re- 
sults suggest and motivate a study aimed at investigating 
the properties of the continuum emission of quasars in the 
UV to NIR region, in order to compare and contrast the 
SEDs of radio loud and radio quiet objects. 

In Section [2] we focus on the quasar sample selection 
criteria and on the host galaxy contribution. In section [3] 
we consider the SED construction and discuss the spectral 
shape. The RLQ and RQQ SEDs are then compared and 
contrasted. In the last section we provide a summary and a 
discussion of the results. 

Throughout this paper, we adopt a concordant cosmol- 
ogy with Ho = 70 km s"^ Mpc"\ rim = 0.3 and Q.k = 0.7. 



2 SAMPLE SELECTION 

In order to mark similarities and differences of the contin- 
uum spectroscopic properties of RLQs and RQQs, our first 
aim is the selection of a quasar sample which is minimally 
biassed against the radio properties and the nuclear color of 
the objects. We start our analysis from the Sloan Digital Sky 
Survey (SDSS) quasar catalogue HI (Schneider et al. 2005); 
we then require that objects are in the observation area of 
the Faint Images of the Radio Sky at Twenty-cm (FIRST, 
Becker et al. 2003), which is the reference catalogue for ra- 
dio data, in order to allow the distinction between RLQs 
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and RQQs. Finally, this sample is cross-correlated with the 
Two Micron All Sky Survey (2MASS, Cutri et al. 2003) in 
order to have detections in the near-infrared (NIR) bands. 



corrected for Galactic extinction, following the indications 
given in Schneider et al. (2005). 



2.1 SDSS quasar catalogue 

The SDSS quasar catalogue (46420 objects) covers about 
4188 deg^ and selects objects that have luminosities larger 
than Mi = —22, have at least one emission line with FWHM 
larger than 1000 km/s or are unambiguously broad absorp- 
tion line objects, are fainter than i = 15.0 and have highly 
reliable redshifts. 

Note that the SDSS quasar catalogue suffers from a 
small bias in favour of radio loud objects, as a FIRST detec- 
tion is one of the starting criteria to select quasar candidates; 
anyway White et al. (2007) verified that this introduces at 
most a very minor bias towards higher radio-optical ratios. 



2.2 FIRST: distinction between RLQs and RQQs 

Starting from the SDSS quasars, we focus on those which 
are in the FIRST observation area, but we don't exclude a 
priori the objects which are under the FIRST flux limit (1 
mjy), so that no bias against the radio power is introduced. 

The distinction between radio loud and radio quiet ob- 
jects is usually made referring to the ratio between the radio 
and optical flux. Here we assume that a QSO is radio loud 
if: 

_ F.(1.5 ■ lO^Hz) 

and radio quiet otherwis^H (e.g. Kellermann et al. 1989). 
Note that our division between RLQs and RQQs should 
not be considered a physical bimodality, an issue on which 
there is still controversy in the literature (e.g. Goldschmidt 
et al. 1999; Ivezic et al. 2002; Jiang at al. 2007). Yet our divi- 
sion is a simple separation in two groups of high or low radio 
power with respect to a certain limit. We verified that all 
the results presented here remain substantially unchanged 
adopting different limits (5 or 30) for the distinction be- 
tween RL and RQ objects. 

There are not radio data for all the objects in the sample 
(91% of the objects are below the FIRST radio fiux limit), 
and obviously the larger part of these are expected to be 
radio quiet objects. In order not to introduce a bias against 
the ratio of RLQs and RQQs, we select all the objects which 
have g < 18.9, so that for this subsample we can discriminate 
between RLQs and RQQs. In fact, if an object has a radio 
flux under the FIRST limit (i.e. F^(1.5 ■ lO^Hz) < ImJy) 
and g < 18.9 (i.e. F^{6 ■ lO^^Hz) > O.lmJy), because of 
condition ([l]) it is a RQQ, while we could not discriminate 
if g > 18.9. This restriction reduces the sample by 66%, 
keeping the original ratio between the number of RQQs and 
RLQs. We obtain a sample of 14395 objects, 1105 of which 
are RLQs and 13290 are RQQs. From the SLOAN data, 
we derived the u, g, r, i, z magnitudes. The values were 



2.3 NIR detections from 2MASS 

The SLOAN-FIRST sample was then cross-correlated with 
the 2MASS catalogue, in order to obtain the J, H and K 
magnitudes of the objects. The 2MASS survey is a collection 
of near-infrared uniformly-calibrated observations of the en- 
tire sky; sources brighter than about 1 mJy in each band 
were detected with a signal-to-noise ratio greater than 10. 
The cross-correlation strongly reduces the sample, which 
now consists of 1761 common objects with observations in 
8 bands each. We assumed two objects are matched if their 
separation angle is less than 0.2 arcsec, which roughly cor- 
responds to the position error of the catalogue^. Obviously 
the request of a 2MASS detection introduces a bias in the 
sample, as among high redshift (low luminosity) objects the 
redder ones are more likely detected. In the following we will 
come back on this issue. 



2.4 Host galcixy component 

A host galaxy component could affect the luminosity 
measures in part of the selected quasars. We require that 
the selected objects to have a negligible galaxy component 
with respect to the nuclear flux. Assuming that all the 
galaxies are ellipticals (SED from Mannucci et al. 2001), 
purely passively evolving (Bressan, Chiosi & Fagotto 1994), 
with Mr = -22.6 for RLQs and Mr = -21.8 for RQQs 
at 2; ~ (see for example Kotilainen et al. 2006), we 
give a flrst order estimate of the galaxy component flux 
in the 8 SLOAN & 2MASS fiffers. An object is excluded 
from the sample if the estimated flux of the host galaxy is 
greater than 20% of the observed flux. The excluded objects 
(~ 1000) are mostly weak, nearby (2 < 0.5) quasars. We 
also tested if the removal of quasars with substantial host 
galaxies could be done by looking at the difference between 
PSF and model magnitudes in SDSS. The results reported 
in the following remain substantially unchanged. 



2.5 Final sample 

The final sample (see Table [T] available in complete form 
electronically) consists of 887 QSOs, of which 774 RQQs 
and 113 RLQs. 

Figures [1] and [2] show the magnitude and redshift distri- 
butions of the radio loud and radio quiet subsamples, indi- 
cating a substantial difference between the two classes. This 
behaviour is well known and it is usually attributed to the 
cosmological evolution of quasar radio properties or to a dif- 
ferent density evolution for RLQs and RQQs (e.g. White et 
al. 2007; Jiang et al. 2007). The issue is obviously relevant 
to our discussion and it is considered in section |3] 



^ Equation U} should be applied to rest-frame data, but this cor- 
rection requires the assumption of a spectral shape of the objects. 
Here we choose to evaluate this ratio based on observed data; we 
verified that the difference is completely negligible. 



^ We verified that this limit on the separation angle causes the 
omission of about 20% objects which belong to both catalogues, 
but drastically reduces the probability of mismatch in the selec- 
tion of the objects. 
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Table 1. The final sample. Monochromatic fluxes are derived from SDSS (col. 4-8), 2MASS (col. 9-11) and FIRST (col. 12) and corrected 
for Galactic extinction; units are 10~^Jy. In col. (12), "n.d." means "not detected". The complete version of this table is available in 
electronic form. 
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Figure 1. Cumulative redshift distribution of the final sample. 
Dotted line refers to the RLQ and dashed line to the RQQ. 



Figure 2. Cumulative Mji distribution of the final sample. Dot- 
ted line refers to the RLQ and dashed line to the RQQ. 



3 SPECTRAL ENERGY DISTRIBUTIONS 
(SEDS) 

In order to construct and compare the SEDs of RQQs and 
RLQs, we first evaluated the rest-frame SEDs of single ob- 
jects. For each one, the SED consists of 8 data-points log(j^)- 
log(i/L^), which are linearly interpolated. In order not to be 
affected by the well known bias which favours more lumi- 
nous objets at higher redshift, using the standard procedure 
(e.g. Elvis et al. 1994), the SEDs of the radio loud and ra- 
dio quiet samples have then been normalized separately at 
10"-«Hz. 

The average SEDs of the two samples have been de- 



termined evaluating the average value of \o%(yLv) at each 
frequency. 

Fig. E] shows mean SED of the whole sample and those 
of the RLQ and RQQ subsamples, and in Table [2] we give 
the number of data-points which contribute to the mean 
for selected values of log(i'), the corresponding average z 
and the mean value of \Q%(yLi,) with corresponding errors. 
The overall SEDs are similar to those reported by Elvis et 
al. 1994 and Richards et al. 2006 (see Fig.gl. The strongest 
feature is the so-called 1 micron inflexion, where the slope of 
the SED changes sign in a vLi, vs. v representation. Here, 
the optical-ultraviolet continuum rises above the infrared 
and forms a ''UV bump" (e.g. Shields 1978; Malkan & Sar- 
gent 1982; Elvis et al. 1994), which is most often interpreted 
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Figure 3. SED for the entire sample of quasar (solid line), RLQs 
(dotted line) and RQQs (dashed line). The horizontal lines show 
the spectral region available at different redshifts. The grey re- 
gions indicate the spectral ranges where the power law fit is per- 
formed. 
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Figure 4. SED for our entire sample of quasar (thick solid line), 
compared to the SEDs obtained by Elvis et al. 1994 (thin solid 
line) and by Richards et al. 2006 (thin dashed line). The three 
SEDs have been normalized at lO^*'^ Hz. 



in terms of thermal emission from an accretion disc (e.g., 
Malkan 1983; Czerny & Elvis 1987). 

At high frequencies (y > lO^'^-^^Hz) the SEDs of all the 
objects are strongly contaminated by the Lyman a forest; 
the small bump at ^ 10^'* ''^ Hz is due to the presence of the 
Ha large emission line. 

It is apparent that RLQs are softer than RQQs. In par- 
ticular, we compare the slope between the 1 micron inflexion 



and the Lyman a forest. The spectral index a (such that 
oc I/") is derived for each object by fitting with a power 
law the observed fiuxes in the spectral range between lO^**'^ 
and 10^^'^^ Hz, excluding the data points affected by the 
most prominent emission features (i.e. ~ 10^^ '' — 10^*'^^ Hz, 

1014.9 _ IQlS l Hz). 

Table [3] reports the mean values of a, the standard de- 
viation and the error on the mean value for the RL and 
RQ subsamples and Figure [5] shows the cumulative distri- 
bution of the spectral index for the two populations. The 
slopes of the RLQs and the RQQs are statistically different 
at more than 99% confidence level, with RLQs being softer. 
The spectral indices of the two populations in the optical re- 
gion differ by ~ 0.2, with a = -0.55±0.41 (±0.04) for RLQs 
and a = -0.31 ± 0.32 (±0.01) for RQQs (where the uncer- 
tainties are the standard deviations and the errors on the 
mean respectively). A Kolmogorov-Smirnov test indicates 
that the probability of RLQs and RQQs to be drawn from 
the same underlying population is negligible [P ~ 10~^). 
This is a strong indication that the two populations are in- 
trinsically different, in the sense that radio loud objects are 
systematically redder than radio quiet quasars. 

In Section 12.31 we noted that the request of a 2MASS 
detection introduces a bias in the sample, as among high 
redshift (low luminosity) objects the redder ones are more 
likely detected. Consequentely, the blue side of the result- 
ing SEDs (which is dominated by high redshift objects) are 
slightly softer with respect to the SED of an unbiassed sam- 
ple. This effect can be quantified by constructing the mean 
SEDs on the whole SLOAN sample, i.e. without the request 
of a 2MASS detection. The shape of the SEDs is modified 
as expected by the bias introduced, but it is noticeable that 
the difference between the SEDs of RLQs and RQQs (which 
is our main result) remains substantially unchanged if we 
consider the whole SLOAN sample to construct the SED 
for log(i^) > 14.8 (where the 2MASS data are uninfiuent) 
and the SLOAN-2MASS sample to construct the SED for 
log(!/) < 14.8 (where the introduced bias is negligible). The 
average spectral indices of the RL and RQ populations, as 
determined on these SEDs, are respectively a = —0.40 and 
a = -0.22, with again Aa ~ 0.2. 

Since we are interested in the study of the shape of the 
SED with reference to the radio power, independently of 
any other parameter, we now construct a number of RQQ 
samples matched in absolute magnitude to the RLQ^. We 
verify that the RQQ matched samples automatically result 
well matched to the RLQs also in redshift. 

Table [3] reports the average results of 10 random 
matched samples. The SEDs of the RQQ matched samples 
are on average slightly softer than those of the entire RQQ 
population (Aa — 0.05), however the Kolmogorov-Smirnov 
test indicates that the probability of RLQs and all the RQQ 
matched samples to be drawn from the same underlying pop- 
ulation is still negligible {P ~ 10~^; the value is greater than 
before as the sample has been drastically reduced), and the 
difference of the mean spectral indices is again ~ 0.2. This is 
consistent with the fact that the whole RQQ population and 



An RQQ with a given luminosity L is extracted from the RQQ 
population with probability proportional to the luminosity den- 
sity of the RLQs. 
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Table 2. For selected frequencies, we report the number of data-points which contribute to the SED of RLQs and RQQs, the mean 
redshift and its standard deviation, the luminosity, its standard deviation and the error on the mean. The luminosities are normalized 
at log{uL^{10'^^-^Hz)). 
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Table 3. Mean values of the spectral index a in the optical— f/V region, standard deviation and error on the mean, for the RL and RQ 
subsamples. Col. 6 reports the Kolmogorov-Smirnov probability that the RLQs and the RQQs (or RQQ matched sample) are drown 
from the same population. 
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an RQQ matched sample have identical distributions in a 
{P ~ 90%) and it is a strong indication that the radio loud 
and radio quiet populations are intrinsically different, in the 
sense that radio loud objects are systematically redder than 
radio quiet quasars independently of the average luminosity 
or cosmic time. 

Finally we compare RLQs and RQQs in different red- 
shift bins (see again Table |3]), and this test shows that the 
color difference between the RL and RQ populations is ap- 
parent in all the redshift ranges. 



4 SUMMARY AND DISCUSSION 

The aim of the present paper is the study of the contin- 
uum emission in the NIR to UV region of a quasar sample, 
in order to compare and contrast the SEDs of the RLQ and 



RQQ subsamples. We selected a sample of quasars with both 
SLOAN and 2MASS detection, to study a spectral range 
from the NIR (for low-z objects) to the UV (for high-z 
QSOs). The sample consists of 887 objects, of which 113 
RLQs and 774 RQQs. 

For each subsample we constructed the mean SED and 
evaluated the average spectral index in the NIR to UV re- 
gion. The slope of the underlying power-laws is significantly 
different: the spectral indices of the two population are sta- 
tistically different at more than 99% confidence, both for the 
whole sample and dividing the sample in redshift bins. The 
difference is present also considering samples well matched 
in absolute magnitude and redshift. 

If the blue bump is due to the superposition of black 
body emissions from an accretion disc, then the color dif- 
ference between RLQs and RQQs should be interpreted in 
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Figure 5. Cumulative distribution of tlic spectral index for the 
RL (dotted line) and RQ (dashed line) subsamples. 
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Figure 6. Comparison of the SED of RQQs (dashed line) to 
the SED of RLQs (dotted line). When RLQs are corrected for 
an additional dust extinction (AA\/ = 0.16mag, solid line), no 
difference is apparent. 



terms of diflerent mean temperatures, in the sense that 
RQQs are hotter. 

We first examine the possibihty that the difference in 
the SEDs of RLQs and RQQs is due to an enhanced dust 
extinction in radio loud objects, as it has been suggested 
e.g. by Francis et aL (2000). In Fig. [6]we compare the SED 
of RQQs to the SED of RLQs, when RLQs are additionaUy 
corrected for dust extinction with respect to RQQs: the off- 
set between the SEDs of the two population is minimized 
adopting ^Ay = O.lGmag. The differences are now com- 



pletely negligible, supporting the hypothesis that RL ob- 
jects are more subject to dust extinction than RQQs. The 
problem would be reconducted to understand why RLQs are 
more extincted than RQQs. A prion this could be related 
to a difference in the inclination angle distribution. Alter- 
natively, it could be that the conditions of dust production 
are related to those justifying large radio emission. 

We then focus on the possibility that the difference of 
disc temperature of RLQs and RQQs is real. Since the tem- 
perature of the inner disc scales as Afg^"'^^ (e.g. Shakura & 
Sunyaev 1973), the difference may be attributed to the fact 
that the black holes of radio loud quasars are supposedly 
more massive (e.g. Dunlop et al. 2003; Falomo et al. 2004; 
Labita et al. 2007). 

The color difference may be linked to the BH spin 
(Stawarz, Sikora & Lasota 2007), as radio emission is usu- 
ally ascribed to a faster spinning. However, spinning BHs 
are expected to have a shorter last stable orbit radius, and 
then a hotter disc, inconsistently with our results. 

Obviously it is possible that the non-thermal (power 
law) continua of RLQs and RQQs are intrinsically differ- 
ent. Supposing that the non-thermal component accounts 
for 80% at lO^ '^Hz, while at lO^'^'^^Hz (in the Big Blue 
Bump) the thermal component is dominant and accounts for 
80%, our data are consistent with a picture where the un- 
derlying power-laws of RLQs and RQQs have qrlq = —1.2 
and aRQQ = — 1.0 respectively, whereas the thermal bumps 
are indistinguishable. The problem would be reconducted 
to understand why the non-thermal continuum of RLQs is 
softer than RQQs. It has been suggested (i.e. Francis et 
al. 2000) that in radio loud samples there is a significant 
chance of synchrotron contamination of the rest-frame R- 
band nuclear luminosities, due to the presence of the rela- 
tivistic jets. This effect explains the red color of QSOs in 
high-frequency selected radio samples (i.e. the PKS survey, 
e.g. Francis et al. 2001), that suffer from a bias towards pole- 
on radio sources which are relativistically boosted above the 
survey flux limit, but no obvious explanation can be invoked 
to justify the color diflerence in our sample. 
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